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Abstract. 

We apply Dynamical Mean-Field Theory to strongly interacting fermions in an 
inhomogeneous environment. With the help of this Real-Space Dynamical Mean-Field 
Theory (R-DMFT) we investigate antiferromagnetic states of repulsively interacting 
fermions with spin ^ in a harmonic potential. Within R-DMFT, antiferromagnetic 
order is found to be stable in spatial regions with total particle density close to one, 
but persists also in parts of the system where the local density significantly deviates 
from half filling. In systems with spin imbalance, we find that antiferromagnetism is 
gradually suppressed and phase separation emerges beyond a critical value of the spin 
imbalance. 
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1. Introduction 

Ultracold atoms in optical lattices provide a new laboratory for interacting quantum 
many body systems Bosonic atoms in optical lattices realize the Bose-Hubbard 
model [2] and undergo a superfluid-Mott insulating transition when the potential 
depth of the optical lattice is increased [3]. Recent experiments directly observed 
correlated particle tunneling [3] and superexchange [5], which are the basic mechanisms 
underlying quantum antiferromagnetism. Also the Fermi surface of fermionic atoms 
in optical lattices |6] and fermionic superfluidity of attractively interacting lattice 
fermions [7j were recently observed, bringing the realization of strongly correlated 
many-fermion states nearby [8]. A two-component mixture of repulsively interacting 
fermions (e.g. ^Li or ^°K) at half filling is predicted to form a correlated paramagnetic 
Mott insulator state above the critical (Neel) temperature and antiferromagnetic order 
below this temperature [9]. Reaching this temperature is predicted to be within today's 
experimental capabilities [TT] . 

In contrast to solid state systems, lattice defects, impurities and phonons are absent 
in optical lattices. However, the spatial inhomogeneity due to the harmonic confinement 
potential is always present, leading to a spatially varying local density. Therefore, 
the concept of long range order is questionable and ordered phases are expected to 
develop on finite length scales. A Hartree-Fock static mean-field theory predicts that 
antiferromagnetism, with staggered magnetization on a finite length scale, coexists with 
paramagnetic states in various spatial patterns, e.g. antiferromagnetism in the center of 
the trap or antiferromagnetism in a ring surrounded by a particle- or a hole-doped atomic 
liquid [12]. On the other hand, both commensurate and incommensurate spin-density- 
waves have been predicted for the hole-doped Hubbard model [131 El IIS] • However, the 
existence and properties of any ordered state on a finite length scale are strongly sensitive 
to quantum and thermal fluctuations. Therefore a theoretical description that captures 
effects of strong correlations and spatial inhomogeneity in a unified framework is needed. 
In this paper we apply a Real-Space Dynamical Mean-Field Theory (R-DMFT), which 
is a comprehensive, thermodynamically consistent and conserving mean-field theory for 
correlated lattice fermions in the presence of an external inhomogeneous potential. The 
R-DMFT takes into account local correlations exactly [16], [HI [IHl [19]. 

We prove that for spin-i lattice fermions with local repulsive interaction 
antiferromagnetic order exists at zero temperature when the harmonic potential is 
present. We find that antiferromagnetic order is stable in spatial regions with total 
particle density close to one, but persists also in parts of the system where the local 
density significantly deviates from half filling. We also show that for strong repulsion 
phase separation occurs in imbalanced mixtures, when the difference in the particle- 
number of the spin components is large. For weaker repulsion a strong imbalance 
destroys the antiferromagnetic order, but does not lead to phase separation. These 
results are especially intriguing with respect to recent experiments on attractively 
interacting fermions with spin imbalance, which have led to still unresolved questions 
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Figure 1. Real-space magnetization profiles for U = 10 on a square (30x30) lattice; 
&) V — 0.1 and ji^ = jii = ^\h) V — 0.2 and = /^^ = 15. Energies are expressed in 
units of the hopping parameter J. 

regarding the nature of the observed phase separation j20l EH [22]. 

2. Model 

Repulsively interacting fermions in an optical lattice almost perfectly implement the 
Hubbard Hamiltonian 

n = -J^ cl^Cja + UY^ n.^n^i + ^{Vi - fi^)ni^, (1) 

(ij),(T i icr 

where rijo- = clo-Cjo-, and Ci^ (cj^) are fermionic annihilation (creation) operators for an 
atom with spin a at site i, J is the hopping amplitude between nearest neighbor sites 
{ij), [/ > is the on-site interaction, fi^ is the (spin- dependent) chemical potential and 
Vi = Vorf is the harmonic confinement potential. Moreover we define p, = |(/it + /^i) and 
Afi = fi^ — fi^. The parameters of this model are tunable in experiments by a change of 
the lattice amplitude and via Feshbach resonances p]. In the following, J = 1 sets the 
energy unit and we take the lattice constant to be a = 1. 

3. Method 

To obtain the ground state properties of this system, we apply a real-space extension 
of dynamical mean-field theory (DMFT) [161 [13 [HI [191 |23l |2l|. Within R-DMFT 
the self-energy is taken to be local, which is exact in the infinite dimensional limit 
[231 [21] • However, in an inhomogeneous system it depends on the lattice site, i.e. 
Hija = Ha^'Sij, where 6ij is a Kronecker delta. Formerly, a similar scheme has been 
developped for systems with inhomogeneity in one direction [25] . Only recently, models 
with full inhomogeneity have been investigated, in particular the Falicov-Kimball model 
[2S1[2Z], disorderded systems [28j and paramagnetic states of cold fermionic atoms [29j. 
In the R-DMFT method, the Hamiltonian is mapped onto a set of single site 
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problems. The physics of lattice site i is described by the local effective action 



S^= - j dr j dr'Y,4^{r)g^\a. 



~U j dTni^{T)n,^ir), (2) 

which explicitly depend on the site index i. Here r is imaginary time. The function 
go\a,T — t') is a local non-interacting propagator interpreted as a dynamical Weiss 
mean-field which simulates the effect of all other sites [18] and is determined self- 
consistently as follows: firstly, given the local self-energies (cr, iujn) obtained from 
solving the action ([2]), the interacting lattice Green's function follows from the Dyson 
equation in real-space representation 

G(o-, iUnY^ = Go(o-, iuJnY^ - S(o-, iUn), (3) 

where a boldface notation indicates that the quantities are matrices labeled by site 
indices i and j and a;„ are the Matsubara frequencies. The non-interacting lattice 
Green's function is given by Go(o", iujn)~^ = {fJ^a + ^i^n)l — J — V, where 1 is the unity 
matrix. The matrix elements Jij are hopping amplitudes for a given lattice structure 
and Vij = 6ijVi represents a spatially varying potential. Secondly, the diagonal elements 
of the lattice Green's function are identified with the interacting local Green's functions, 
i.e. G^*) {(7,iuj„) = Gii{a,iiUn)- Finally, the Weiss mean-field is obtained from the local 
Dyson equation 

which closes the set of the self-consistency equations. 

The most difficult step in this procedure is the solution of the local action ([2]). 
This step is, however, similar to the solution of the local action in a homogeneous 
DMFT calculation. The difference is that in the present case the Weiss field Qq^ {a, r) 
is obtained via the Real-Space Dyson equation ([3]), which incorparates the effect of the 
spatial inhomogeneity. This implies that for the numerical solution of the local action 
we can use standard techniques, which have proven to be reliable and efficient. In the 
present manuscript we use the Numerical Renormalization Group (NRG) at T = 
[33| [Ml [35| [361 [37] to solve the single site problems. 

In practice the self-consistent solution is obtained iteratively from the initial Weiss 
mean-fields QQ\a,iuJn) which are chosen differently for different spin a and lattice sites 
i. Then the solutions with staggered magnetization or phase separation are obtained 
naturally in contrast to the standard DMFT, where additional sublattice structure has 
to be added p!8]. 

Within R-DMFT significantly larger systems can be investigated than those studied 
by quantum Monte Carlo [30l [3T| [32] . for which in two and three dimensions only 
homogeneous data are available. The computational effort scales polynomially with the 
number of lattice sites within R-DMFT. The application of the real-space Dyson 
equation requires a sparse matrix inversion for each frequency, which scales as N^^"^. 
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Figure 2. Particle and spin density profiles determined within the exact solution of 
R-DMFT and within the Thomas- Fermi approximation (TEA) to R-DMFT. The main 
panel shows the local total density rii and the local magnetization rrii as a function of 
a distance from the trap center. The inset shows Ui and along the 2/ = 5 line. The 
parameters in the a) and b) panel are the same as in the respective panels of Fig. [1] 



The number of NRG calculations per R-DMFT-run is linear in A^, but can be kept 
small due to symmetries. Moreover, the solution of the real-space Dyson equation can 
be parallelized over the frequencies and the NRG-calculations can be parallelized over 
the lattice sites. 

4. Results 

We apply this method to spin-| fermions in a two-dimensional square lattice with 
harmonic confinement. In the context of cold atoms a two dimensional system can 
be realized by applying a highly anisotropic optical lattice, which divides the system 
into two-dimensional slices. Although not exact, R-DMFT is expected to be a 
good approximation for the two-dimensional situation at zero temperature, since the 
derivation of the DMFT equations is controlled by the small parameter 1/z = 1/4 on 
the square lattice. 
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Figure 3. Total density and staggered magnetization as a function of effective 
chemical potential obtained within the Thomas-Fermi approximation to R-DMFT for 
two dimensional square (main panel) and three dimensional cubic lattices (inset). Main 
panel: J7 = 10 (diamonds) and U = 20 (circles); inset: U = 30. 



4.1. Balanced mixture 

First we consider the case of an equal mixture of spin- up and down atoms: = ^i- 
We find that antiferromagnetic order is stable in the presence of the inhomogeneous 
harmonic potential. In Fig. [1] we present examples for the spatial dependence of 
the magnetization at different strengths of the confining potential and the chemical 
potentials. In the case that the lattice at the center of the trap is half-filled, 
antiferromagnetism appears in the center of the system (Fig. [1^). When the particle 
density in the center of the trap is higher, antiferromagnetic order forms in a ring 
enclosing a paramagnetic region (Fig. Wp)- These results are particularly important 
for ongoing attempts to realize antiferromagnetic states in optical lattices. Namely, we 
predict that the observation of antiferromagnetic order does not critically depend on the 
number of atoms in the system. For sufficiently strong repulsion between the particles, 
the necessary condition to find antiferromagnetic order is to prepare the system such 
that the local filling factor approximates or exceeds one in at least part of the system. 
We find no evidence for phase separation or a paramagnetic insulating boundary layer 
for the = case. 

The antiferromagnetic ground state of homogeneous fermions described by the 
Hubbard Hamiltonian ([T]) without trap is stable when the density of particles varies 
from n ^ 0.8 to 1.2, depending on the interaction value U [38]. On the contrary, in the 
presence of the external harmonic potential, antiferromagnetic order appears for much 
lower or higher local total densities. Indeed, in Fig. [2] we present examples of the local 
density and the local magnetization rui = ^{uii — n^j) as a function of distance from 
the center (main panel) and along a cut through the system (inset) which proves that 
antiferromagnetic order extends from the center of the trap and disappears only when 
Hi ^ 0.5 in Fig. [2^). Similarly, Fig. [2)d) shows that antiferromagnetic order is stable on 
a ring when the local density extends between 0.5 < < 1.5. 
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We also determine the local density and the local magnetization within the Thomas- 
Fermi approximation (TFA) to R-DMFT, where the external potential is only included 
by a spatially varying chemical potential [39] . The agreement between the full R-DMFT 
and TFA results is very good in regions well within or outside the antiferromagnetic 
domain. Encouraged by this, Fig. [3] shows additional TFA+R-DMFT profiles that 
can be used to compare R-DMFT with experiments for realistic systems in two and 
three dimensions. However, the staggered magnetization decays abruptly within TFA 
as compared to the full R-DMFT solution, i.e. the Thomas-Fermi approximation to R- 
DMFT essentially reproduces results from the standard homogeneous DMFT, cf. Fig. [31 
The wider stability regime of the antiferromagnetic order found within full R-DMFT is 
caused by a proximity effect; antiferromagnetic order is induced into parts of the systems 
where the local densities are too low to stabilize antiferromagnetism in the homogeneous 
case. On the other hand, the proximity of the paramagnetic state reduces the staggered 
magnetization when the local density is close to one. 

4-2. Imbalanced mixture 

We now proceed by investigating the imbalanced case [IQ], i.e. N-^ ^ N^. Imbalance 
between the two spin-components is induced by a nonzero chemical potential difference 
Afi = /i| — which corresponds to a magnetic field. In the experimental situation, 
the density imbalance can be highly tuned and is stable due to the suppression of spin- 
flip scattering proccesses in cold-atomic gases. Representative results are presented in 
Fig. [H where we plot the up- and down-component of the density along a cut through 
the system. Upon increasing the imbalance parameter A/i, we find suppression of the 
antiferromagnetic order and emergence of phase separation between the minority and 
majority species. The phase separation region starts to develop far away from the 
center of the trap at small A/x and gradually spreads toward the center. We thus find 
that the border of the anti-ferromagnetic domain is most sensitive to phase separation. 
This is indeed reasonable: the energy cost to polarize the antiferromagnetic state is the 
energy difference between an antiferromagnetic state and a ferromagnetic state. This 
is of the order J'^/U, which small for the large interaction U considered here. The 
antiferromagnetic order is thereby more unstable for larger distances to the trap-center, 
because of the vicinity to the paramagnetic region. The energy cost to polarize the 
paramagnetic regime is higher, because in this case kinetic energy has to be paid, whereas 
in the anti-ferromagnetic domain the kinetic energy is already quenched because the 
particles are almost localized. Due to the proximity effect we find that the paramagnetic 
regime close to the insulating domain also gets phase-separated, which leads to a ring- 
like structure of the minority species. 

At strong interaction, U = 10 in the case shown in Fig. [H atoms with different spins 
ultimately tend to occupy different spatial regions to avoid the mutual interaction and 
the minority species is completely expelled from the trap center. At weaker interaction, 
however, we found that the imbalanced system still contains interpenetrating atoms 
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Figure 4. Spin resolved particle densities for an imbalanced mixture obtained within 
R-DMFT for U = 10. Panels (a-c) show component densities along the y = \ line for 
gradually increasing imbalanced A/i — 0.3 (a), 0.75 (b), 1 (c). The two lower panels 
show the space resolved up- (d) and down- (e) density for A/i = 1. The lattice size is 
20 X 20 and other parameters are: V — 0.2, p, — 5. 
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Figure 5. Spin resolved particle densities for an imbalanced mixture obtained within 
R-DMFT for U — 7.5. Panels (a-c) show component densities along the y — line for 
gradually increasing imbalanced A/i = 0.4 (a), 0.6 (b), 8 (c). The lattice size is 20 x 20 
and other parameters are: V = 0.2, p, — 3.75. 

with different spins and pliase separation does not occur. Tfiis is sfiown in Fig. wfiere 
for A/i = 0.8 tlie antiferromagnetic order has completely disappeared, but the two spin 
components are still interpenetrating. The small oscillations in the component densities 
can be understood as Friedel oscillation due to the small size of the system. We note 
that in the case of imbalanced spin-mixtures the agreement between the TEA and the 
exact solution to R-DMFT is far less good than in the balanced case presented above. 

5. Conclusions 

In conclusion we used the R-DMFT to establish stability of antiferromagnetism for 
balanced fermionic spin-| systems in a trap and the appearance of phase separation 
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for imbalanced mixtures. The antiferromagnetic order predicted here can be observed 
at low enough temperatures by Fourier-samphng of time-of-flight images via Raman 
pulses [H], by measuring spin correlation functions via local probes [12], probing noise 
correlations [ISj Hi], polarization spectroscopy [15], and Bragg scattering [16]. The 
effect of spatial inhomogeneity on these probes will be investigated within R-DMFT in 
future studies. Moreover, the R-DMFT scheme presented here opens up the possibility 
to study a variety of other strongly correlated systems in inhomogeneous environments. 
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